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Stellar Evolution with Rotation and Magnetic Fields 
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Abstract. We examine the generation of a magnetic field in a solar-like star and its effects on the internal dis- 
tribution of the angular velocity. We suggest that the evolution of a rotating star with magnetic fields leads to 
an equilibrium value of the differential rotation. This equilibrium is determined by the magnetic coupling, which 
favours a constant rotation profile, and meridional circulation which tends to build differential rotation. The 
global equilibrium stage is close to solid body rotation between about 0.7 and 0.2 i?0, in good agreement with 
helioseismic measurements. 
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1. Introduction 

One of the most severe problems in stellar physics con- 
cerns the rotation profile of the radiative interior of the 
Sun. Helioseismological results indicate that the angular 
velocity i^{r) is constant as a function of the radius r be- 
tween about 20 % and 70 % of the total solar radius Rq 
(Brown et al. HM^ Kosovichev et al. [TM7I Couvidat et 
al. I2003|l . while meridional and rotational turbulent dif- 
fusion (e.g. Zahn 1992j produce an insufficient internal 
coupling to ensure solid body rotation (Pinsonneault et 
al. 119891 Chaboyer et al. I1995|l . This suggests that an- 
other effect intervenes. Mestel & Weiss H1987f) proposed 
that a weak internal magnetic field could provide the 
required internal coupling. Charbonneau & MacGregor 
H1993[l also showed that large-scale internal magnetic fields 
can yield a weak internal differential rotation at the so- 
lar age. Another proposition to explain the flat rotation 
profile of the Sun is the angular momentum transport by 
internal gravity waves (Zahn et al. 119971 Talon et al. 120021 
Talon & Charbonnel lMTSji . 

An efficient dynamo has been proposed to operate 
in stellar radiative layers in differential rotation (Spruit 
I2002|l . This dynamo is based on the Tayler instability, 
which is the first one to occur in a radiative zone (Tayler 
[1^751 Pitts & Tavler ligHBj) . Even a very weak horizontal 
magnetic field is subject to Tayler instability, which then 
creates a vertical field component, which is wound up by 
differential rotation. As a result, the field lines become 



progressively closer and denser and thus a strong horizon- 
tal field is created at the energy expense of differential 
rotation. 
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Maeder & Meynet studied the effects of the 

Spruit dynamo on the evolution of massive stars. They 
showed that a magnetic field can be created during the 
main sequence evolution of a rotating star by the Spruit 
dynamo and examined the timescale for the field creation, 
its amplitude and the related diffusion coefficients. Maeder 
& Meynet l|2004|l also developed a generalisation of the dy- 
namo equations in order to encompass all cases of fi- and 
T-gradients, as well as all cases from the fully adiabatic to 
non-adiabatic solutions. The clear result of these studies 
is that magnetic field and its effects are quite important 
in massive stars. The magnetic coupling, resulting from 
the Tayler-Spruit dynamo, is found to be able to enforce 
solid body rotation in massive stars. The question is now 
whether this dynamo also operates efficiently in a slow ro- 
tator like the Sun, so as to be responsible for the observed 
constancy of ri(r). 

In Sect. 2, we collect in a short consistent way the basic 
equations of the dynamo. In Sect. 3, the numerical models 
are presented, while Sect. 4 gives the conclusion. 



2. The dynamo equations 

In this section we briefiy summarize the consistent system 
of equations for the dynamo (see Spruit .2002:. and Maeder 
& Meynet Unni for more details). 
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The energy density ub of a magnetic field of intensity 
B per volume unity is 

B 



B^ 1 2 2 

— = - pr ujj^ with UJA 

O TT Z 



(4 TTp) 2 r 

where uja is the Alfven frequency in a spherical geometry. 
In stable radiative layers, there is in principle no partic- 
ular motions. However, if due to magnetic field or rota- 
tion, some unstable displacements of vertical amplitude 
1/2 occur around an average stable position, the restoring 
buoyancy force produces vertical oscillations around the 
equilibrium position with a frequency equal to the Brunt- 
Vaisala frequency N (see Eq. 14 of Maeder & Meynetl2U041 
for the definition of iV in a medium with both thermal and 
magnetic diffusivity K and rj). 

The restoring oscillations will have an average density 
of kinetic energy 



(2) 



where /n is a geometrical factor of the order of unity. If 
the magnetic field produces some instability with a ver- 
tical component, one must have ub > wn- Otherwise, 
the restoring force of gravity which acts at the dynamical 
timescale would immediately counteract the magnetic in- 
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stability. From this inequality, one obtains P < j^- 

If /n = 5 , we have the condition for the vertical amplitude 
of the instability fSoruit lMTH Eq. 6), 



I < 



N 



(3) 



where r is the radius. This means that there is a maximum 
size of the vertical length I of a magnetic instability. In 
order to not be quickly damped by magnetic diffusivity, 
the vertical length scale of the instability must satisfy 



(4) 



where fl is the angular velocity and as the characteris- 
tic growth-rate of the magnetic field. In a rotating star, 
this growth-rate is ctb — (^a/^) due to the Coriolis force 
(Spruit 2002! see also Pitts & Tavler [T?M|l . The combina- 
tion of the limits given by Eqs. Q and Q gives for the 
case of marginal stability. 
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rj 
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(5) 



The equality of the amplification time of Tayler insta- 
bility Ta = N/{ujA^q) with the characteristic frequency ctb 
of the magnetic field leads to the equation C Spruit l^002|l 

T^^^TT With q = -^, ■ 

iZ N oinr 

By eliminating the expression of iV^ between Eqs. (jSJ 
and ©, we obtain an expression for the magnetic diffu- 
sivity. 



(6) 



(7) 



Eqs. (0) and © form a coupled system relating the two 
unknown quantities 77 and loa- Instead, one may also con- 
sider for example the system formed by Eqs. © and (0). 
Formally, if one accounts for the complete expressions of 
the thermal gradient V, the system of equations would be 
of degree 10 in the unknown quantity x = (^f ) (Maeder 
& Mevnet l[2UIMI . The fact that the ratio -q/K is very smah 
allows us to bring these coupled equations to a system of 
degree 4 (Maeder & Meynet [ 



(iV| + Nl) - 



K 



x^ + 2N^x- 2n^q^ = . (8) 



The 



K is the radiative diffusivity, namely K — ^^'i^^j 
solution of this equation, which is easily obtained numer- 
ically, provides the Alfven frequency and by Eq. Q the 
thermal diffusivity. 

The azimuthal component of the magnetic field is 
much stronger that the radial one in the Tayler-Spruit 
dynamo. We have for these components f Spruit f5()()2|) 



B^ = (47rp) 2 r uja and B^ = B^ (IJr) , 



(9) 



and 



where uja is the solution of the general equation 
is given by Eq. |j2Jl. 
Turning towards the transport of angular momentum 
by magnetic field, we first write the azimuthal stress by 
volume unity due to the magnetic field 



S — - — ByB^p — - — 

4 TT ^ 4 TT 



Bl^ pr^ 



N 



(10) 



Then, the viscosity 1/ for the vertical transport of angular 
momentum can be expressed in terms of S (Spruit i20C 



S 



p q fl 



\n ) [n 



(11) 



This is the general expression of v with loa given by the 
solution of Eq. . We have the full set of expressions nec- 
essary to obtain the Alfven frequency lua and the magnetic 
diffusivity rj. Let us recall that rj also expresses the verti- 
cal transport of the chemical elements, while the viscosity 
v determines the vertical transport of angular momentum 
by the magnetic field. 

3. Stellar models 

We consider here models of 1 Mq , with the solar chemical 
composition of Grevesse & Sauval (1998 ). The stellar evo- 
lution code used for these computations is the Geneva code 
including shellular rotation (Meynet & Maeder r2005|l . We 
use the braking law of Kawaler H1988|l in order to repro- 
duce the magnetic braking undergone by low mass stars 
when arriving on the main sequence. Two parameters en- 
ter this braking law: the saturation velocity Jlgat and the 
braking constant K. Following Bouvier et al. H1997fl . llgat 
is fixed to 14 Hq and the braking constant K is calibrated 
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Fig. 1. Rotation profiles as a function of the lagrangian 
mass in solar units for models with (continuous lines) and 
without magnetic field (dashed lines). The dotted line in- 
dicates the initial solid body rotation on the ZAMS. The 
other lines correspond to an age of respectively 1, 2, 3, 4 
and 4.57 Gyr. fl increases at the centre during the evolu- 
tion on the main-sequence. No magnetic braking at the 
surface is included. 

on the Sun. It is beyond the scope of the present paper to 
compute complete solar models reproducing the detailed 
solar structure deduced from helioseismic measurements. 
In this study, we simply focus on the evolution of the ro- 
tation profile of IMq models which reproduce the solar 
luminosity and radius at the age of the Sun (4.57 Gyr). 

To investigate the effects of magnetic fields on the ro- 
tation profile of a slow rotating solar-like star, we first 
consider models without magnetic braking at the surface. 
Two models are computed: one with rotation only and a 
second with both rotation and magnetic field. The initial 
velocity of these models is approximately equal to the solar 
surface rotational velocity. Fig. ^ shows the evolution of 
the internal rotation profile for both models, starting from 
solid body rotation on the ZAMS. The rotation profile of 
the model with rotation only changes with time during 
the main sequence evolution due to the transport of an- 
gular momentum by circulation. As a result, the model 
with only rotation shows after 4.57 Gyr an angular veloc- 
ity f2 which is monotically increasing when the distance to 
the centre decreases. The situation is quite different when 
magnetic fields are accounted for. As shown in Fig. ^ the 
angular velocity Q is almost constant between the surface 
and about 0.3 Mq. In the central parts, the angular veloc- 
ity increases due to the decrease of the horizontal coupling 
insured by the magnetic field strength resulting from the 
jU-gradient in the stellar core varies like V~^). 



Fig. 2. Same as Fig. but with an initial velocity of 
20kms~^. The dotted line indicates the initial solid body 
rotation on the ZAMS. The other lines correspond to an 
age of respectively 0.1, 0.5, 1, 2, 3 and 4.57 Gyr. The sur- 
face angular velocity decreases when the star evolves due 
to the magnetic braking at the surface 

As a second step, models with initial velocities of 20 
and 50kms~^ are computed. The braking law of Kawaler 
fl988) is used for these models in order to reproduce the 
solar surface rotational velocity. Fig. El compares the evo- 
lution of the internal rotation profile of models with only 
rotation and with both rotation and magnetic fields, start- 
ing with Vini = 20kms~^ on the ZAMS. We notice that 
the surface rotation velocity rapidly decreases due to the 
magnetic braking at the surface. For the model including 
only rotational effects, this results in a large differential 
rotation reaching a factor of about 20 between the angu- 
lar velocity at the surface and in the stellar core at the age 
of the Sun, in good agreement with the previous results of 
Chaboyer et al. fljlSsj) . but in contradiction with the flat 
rotation profile of the Sun. Fig. [3 shows that the model 
with both rotation and magnetic fields displays an almost 
constant angular velocity throughout the radiative inte- 
rior, with only a small increase of f2 in the central parts 
(for Mr < 0.2 Mq). Fig. 131 better compares the theoret- 
ical rotation profiles of models computed with an initial 
velocity of 50kms~^ to the one deduced from helioseismic 
measurements (Couvidat et al. l?D03p . Note that the rota- 
tion profiles of models with an initial velocity of 50 km 
are very similar to those with = 20kms~^. The only 
difference is that models with Vmi = 50kms^^ exhibit 
slightly faster rotating cores than those computed with 
Willi = 20kms^"'^. Fig. 121 clearly shows that the rotation 
profile of models including both rotation and magnetic 
fields is in good agreement with the helioseismic measure- 
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Fig. 3. Rotation profile for a model with rotation only 
(dashed line) and with both rotation and magnetic field 
(dotted line) at the age of the Sun. The initial velocity is 
50kms^^. The points with their respective error bars cor- 
respond to the angular velocities in the solar radiative zone 
deduced from GOLF+MDI and LOWL data (Couvidat et 
al.lMll. 




BO 
O 



-10 



-15 



f. ' i 

i! i 



_L 



_L 



0.2 



0.4 0.6 



0.8 



Fig. 4. Diffusion coefficients in the model with both rota- 
tion and magnetic fields and v-mi — 50kms^^ at the age 
of the Sun. 



ments, while models with only rotation predict a too fast 
increase of when the distance to the centre decreases. 

Finally, the values of the various diffusion coefficients 
for the model with magnetic fields and Vini — 50kms^^ 
are shown in Fig. 0] The largest diffusion coefficient is v 
which acts for the vertical transport of angular momen- 
tum. The large value of v imposes the nearly constant 
fl in the interior and confirms the dominant role of the 
magnetic field for the transport of angular momentum. 
The numerical value for the azimuthal component of the 
field is of the order of a few 10^ G. We also notice the 
small value of the coefficient of the shear turbulent mixing 
-Dshoar- Indeed, it is of the same order of magnitude as the 



coefficient DcS, which applies to the transport of chemi- 
cal elements by meridional circulation, while in a rotating 
star without magnetic field -Dshoar is generally much larger 
(about 4 orders of magnitude). This small value is a conse- 
quence of the near solid body rotation of magnetic models 
and suggests that, for slow rotating solar-like stars, the ro- 
tation induced mixing is less efficient in magnetic models 
than in models with rotation only. This seems to be in 
agreement with observations. A detailed study of the evo- 
lution of trace elements like lithium and '^He with and 
without magnetic field will be needed to really investigate 
this point. Fig.0]also shows that the value rj/K is always 
very small, which justifies the simplifications made in de- 
riving Eq. ijHJ). 

4. Conclusion 

The main result of this study is that the Tayler-Spruit 
dynamo can account for the flat rotation profile of the 
Sun as deduced from helioseismic measurements. There 
remains however some doubts whether this dynamo is re- 
ally active in stellar interiors, since 3D simulations have 
not yet confirmed the existence and efficiency of this par- 
ticular instability; these simulations show a delicate bal- 
ance between the generation of the instabilities and their 
relaxation to stable configurations (Braithwaite & Spruit 
l2(J(J4f) . It is also worthwhile to recall that magnetic field 
is not the only explanation, since purely hydrodynamical 
stellar models including the transport by internal gravity 
waves constitute another promising alternative. 
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